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ABSTRACT The standard treatment of relativistic thermodynamics is not convenient for 
a systematic treatment of mixtures. It is proposed that a formulation of thermodynamics 
as an action principle may be a suitable approach to adopt for a new investigation of that 
and other problems in thermodynamics and astrophysics. 

This third paper of the series applies the action principle to a study of mixtures of 
ideal gases. The action for a mixture of ideal gases is, in the first approximation, the sum 
of the actions for the components. The entropy, in the absence of gravity, is determined 
by the Gibbs-Dalton hypothesis. Chemical reactions such as hydrogen dissociation are 
studied, with results that include the Saha equation and that are more complete than 
traditional treatments, especially so when gravitational effects are included. A mixture 
of two ideal gases is a system with two degrees of freedom and should exhibit two kinds 
of sound. This is not supported experimentally and the Gibbs-Dalton lagrangian has to 
be modified accordingly. In the presence of gravity the Gibbs-Dalton hypothesis must be 
further modified to get results that agree with observation. The overall conclusion is that 
experimental results serve to pin down the lagrangian in a very efficient manner. It leads 
to a convenient theoretical framework in which many dynamical problems can be studied, 
and in which the incorporation of General Relativity is straightforward. 
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I. Introduction 

Most physical theories that have reached a high level of development are formulated 
as action principles. This is true of General Relativity, but it is not true of astrophysical 
applications of General Relativity. The reason for this historical anomaly has been the need 
to incorporate thermodynamics, a highly developed theory that represents the exceptions: 
thermodynamics is not formulated as a dynamical action principle. 

This is odd, since Gibbs' famous contributions to thermodynamics are entirely based 
on axioms of minimum energy and maximum entropy, stopping just short of introducing 
a lagrangian. That Gibbs' "energy" can be regarded as a dynamical hamiltonian, at least 
in certain contexts, is well known; see Fetter and Walecka (1980). 

That astrophysical problems, within a General Relativistic setting, are subject to a 
treatment based on an action principle was conjectured independently by Bardeen (1970), 
Schutz (1970) and Taub (1954). A concrete example, related to the nonrelativistic theory 
of Fetter and Walecka, was presented by us (Fronsdal 2007). But the absence of an 
established formulation of thermodynamics, suitable for the extension to the domain of 
General Relativity, may have led to a justifiable skepticism about the viability of this 
approach to astrophysical problems. 

Whence the present series of papers, intended to present a view of thermodynamics that 
can be extended to the context of General Relativity, much in the manner in which it was 
done in the last cited paper, but with the ther mo dynamical interpretation more developed 
and better understood. 

Our approach to thermodynamics is characterized by a complete dynamical framework 
based on an action principle. The action principle that desribes adiabatic configurations of 
an ideal gas was known. Here we expand the context to include heat transfer, dissipation, 
mixtures and chemical reactions. 

This paper is preparation for taking up one more of the issues that, it seems to us, 
cloud the traditional approach to relativistic thermodynamics: the treatment of mixtures. 
Mixtures have more degrees of freedom than those envisaged by Tolman's relativistic equa- 
tions (Tolman 1934). We begin by determining the lagrangian for a mixturte of ideal gases. 
This requires experimental input and the first observation that is invoked is the validity 
of the Gibbs-Dalton hypothesis concerning the entropy of mixtures. When chemical reac- 
tions are included we obtain an expression for the degree of dissociation that is formally 
identical to the Saha equation. Taking the latter as a convenient, approximate summary of 
experimental results, we find that it places a very simple constraint on the entropy. Once 
the lagrangian is determined all applications are straightforward. The determination of 
the lagrangian by means of observation is a main strategy adopted in this paper. 

A principal feature of an action principle is that all properties of a system are strongly 
interrelated. Having arrived at a lagrangian that is adequate for some applications, one 
must try to use the same lagrangian in other situations. And if other considerations 
require a modification of the lagrangian, then one has to accept that applications already 
thought to have been done will have to be revisited. For example, having found a way to 
account for some of the observed features of sound propagation in mixtures, we may find 
that difficulties arise when the theory is applied to the problem of concentrations in the 
atmosphere. Any modification of the lagrangian will have effects on all applications. For 
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another example, a polytropic model of the atmosphere cannot remain polytropic when 
the energy of the photon gas is taken into account, as it must for very high temperatures. 

Summary 

Section II is a brief account of the first paper in this series. The principal dynamical 
variables for describing a simple gas are the density, a velocity field and the temperature. 
Variation of the action with respect to the temperature gives a fundamental relation of the 
ideal gas, an expression for the entropy in terms of the density and the temperature. Any 
discussion of mixtures is intimately concerned with entropy. A study of simple processes, 
including cooling and free expansion, reveals just enough of this concept to allow us to 
understand mixtures, including chemical reactions. 

To deal with mixtures we need additional input from experiment or from kinetic theory. 
A key statement that is taken from observation is the Gibbs-Dalton hypothesis about the 
entropy of mixtures. Traditional applications make use of general properties of the Gibbs 
function but all these properties are relative to the Gibbs function of a certain "reference 
configuration" and the investigation is hampered by the fact that this latter function is 
not known. For this reason, the theoretical treatment of reaction rates, for example, has 
limited predictive power and refers instead to a "reaction constant" that is actually a 
function of the temperature, to be determined by experiment. In Section III we use the 
experimental information that is expressed by the Gibbs-Dalton hypothesis to determine 
the lagrangian, thus preparing the way for many types of applications. 

Turning to chemical reactions, we find that other experimental input is required. The 
first case studied, in Section III. 3, is the simplest one, the atomization of a warm hydrogen 
gas. The Euler-Lagrange equations of motion lead directly to a formula for the reaction 
rate in terms of a single entropy parameter. This equation is formally identical to the 
simplest version of the Saha equation (Eggert 1919, Saha 1920). Taking this equation 
as a convenient summary of experimental results, we observe that this parameter is con- 
stant through the region of dissociation. That is, observation makes an especially simple 
statement about the lagrangian. 

One could now, for example, calculate the characteristics of sound propagation in partly 
dissociated hydrogen. Other results that follows immediately from the equations of mo- 
tion are formulas for the internal energy and the entropy. A lagrangian for more general 
chemical reactions is conjectured, Section III. 4. 

In Section IV we use this lagrangian to study the propagation of sound in a mixture. 
One obtains a system of two linked wave equations and two distinct normal modes, with 
different propagation velocities. As far as we know, this phenomenon has been recognized 
in the context of superfluid helium only (Putterman 1974). But this prediction of dual 
sound speeds is contradicted by experiment, which makes it necessary to modify the la- 
grangian that was built on the Gibbs Dalton hypotheses, without necessarily giving up 
that hypothesis, since it applies to equilibria only. 

In Section V we introduce the gravitational field in the usual way, by including the 
gravitational potential in the hamiltonian, in each of the contexts considered and for 
which the lagrangian has been found. Application of the Gibbs-Dalton hypothesis to this 
situation is ambiguous and unsuccessful. We again turn to measurements (in the earthly 
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atmosphere) to determine the entropy of a mixture in the presence of gravity, obtaining 
density and temperature profiles in good agreement with observation. 

In another paper in this series, we shall return to the context that stimulated this 
research, the application of relativistic thermodynamics to stellar structure and develop- 
ment. The generalization of the simplest lagrangian to include the effect of a generally 
relativistic, dynamical metric is straightforward (Fronsdal 2007). 

As in all the papers of this series, viscosity is neglected and the motion is irrotational. 
This limitation is not as strong as it may appear to be; as may be seen from the existence 
of a hamiltonian formulation of a theory of turbulence (Onsager 1949). 



II. The ideal gas* 



II. 1. The extremal principles of Gibbs 

The basis for a close integration of the action principle with thermodynamics is the 
identification of the hamiltonian with the function 



H = F(V, T) + ST + PV, 



(2.1) 



a function of 4 independent variables V, T, S and P: volume, temperature, entropy and 
pressure. The function F is the free energy. 

The states of thermodynamic equilibrium are the points in this 4-dimensional manifold 
at which H is extremal with respect to variations of T and V, 



dH 



dT 



dH 



OF 



dV 



v,s,p dT 



OF: 

t,s,p dV ' 



+ S = 0, 



+ P = 0. 



(2.2) 



(2.3) 



The system is defined by the choice of the free energy function F. When this function has 
been specified the two relations define a 2-dimensional surface in a 4-dimensional manifold 
with coordinates V, T, S, P. 
In the case of an ideal gas, 



and these relations give 
and 



F(V,T) = -TZT\n(VT n ), 

S/n = n + ln(VT n ), 
PV = 1ZT, 



(2.4) 

(2.5) 
(2.6) 



* The development presented here can easily be expanded to apply, for example, to a gas 
with an arbitrary, cubic equation of state (Fronsdal 2010). 
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respectively. The on shell internal energy density is 



U= (F(V,T) + ST) 



= nIZT. 



on shell 



(2.7) 



"On shell" means that the function is restricted to the surface defined by (2.2). The last 
three relations define the ideal gas; we take that as sufficient proof that (2.4) is a correct 
and useful expression for the free energy, and that the expression (2.1) is a promising 
candidate to serve as hamiltonian. 

The free energy is the only potential that is defined from the beginning as a function of 
two natural variables. The internal energy is a "variable" . It becomes a function of V, T 
and S by identification with F + TS, restricted to the V, T, S surface. It is regarded as 
a function of V and S, defined by the Legendre transformation U(V, S) = F + ST and 
elimination of T by use of the relation (2.2). The function F occupies a special position. 

There is a symmetry between energy and entropy, strongly emphasized by Gibbs (1878). 
One can regard the entropy as a function of V, T, H and P, 
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(H-F(V,T) -PV). 



(2.8) 



The conditions that this function be stationary with respect to variation of the variables 
T and V, with H fixed, are exactly the same as (2.2-3), 



dS 
dT 



V,U,P 



~(s+™ 

T\ b+ dT 



v 







(2.9) 



and 



dS_ 
8V 



T,U,P 



_ 1 ,dF 
~ ~T\dV 



+ p) = 0. 



(2.10) 



Gibbs claims that the two statements, 5H = and SS = 0, are equivalent. But that is not 
true in the local extension of the theory. 



2. Local relations 

The local extrapolation of thermodynamics seeks to promote these relations to field 
equations that describe local but stationary configurations. The functions F, H, S and V 
are given new interpretations as specific densities and specific volume. The mass density 
is p = 1/V and densities /, h, s are defined by 

f(p,T) = pF(V,T), h = pH, s = pS. (2.11) 



The Hamiltonian is, 

J d 3 xh, h = pv 2 /2 + P (j) + f(p,T) + sT + P. (2.12) 
The effect of an external gravitational field is represented by the term pcf). 
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Variation with respect to T, with p, S, T treated as independent variables, leads to 

+ s = 0. (2.13) 



dT 



p 

Variation with respect to p, with the mass and the volume fixed, * gives 



h + p d(f/p) 



P ' dp 

The local thermodynamic pressure is defined by 

Af/P) 



= A. (A = constant.) (2.14) 



r dp 

then the last relation (2.14) reduces to 



- f = P, 

T 



grad ^±^ = 0, (2.15) 



In the case of the ideal gas 



f(p,T) = 1lpTln^, p = KpT. 

The on shell value of h is the internal energy density, 

h\on shell = pv 2 /2 + pcj) + riRpT 
and the variational equations reduce to 

S = n\n(VT n ), and grad(V 2 /2 + </> + (n + 1)^t) = 0. (2.16) 

At equilibrium (v = 0), in the absence of external forces (0 = 0), the temperature is 
uniform. 

Remark. The Bernoulli equation, in the stationary case (all fields time independent), is 

pgrad(-f7 2 /2) + pgradcp + grad^pT = 0. 

It should be noted that it is equivalent to (2.16) only if the gas is isentropic, when p n /T 
is uniform (Fronsdal 2010). 



* A complimentary variation of the density, with the mass fixed but the volume not, gives 
the result that the internal pressure p must agree with the external pressure P on the 
boundary. See Section II. 6. Any sub volume is in equilibrium with the pressure excerted 
by the gas outside its boundary; hence P is identified with the thermodynamic pressure p. 



6 



Reversing the roles of h and s we extremize the total entropy 

j &xa> s =^(h-pv 2 /2-f{p,T)-P). 

Variation with respect to T, with p,T and H = h/p treated as independent variables, 
gives the same result, df/dT + s = 0, but variation with respect to p gives an additional 
constraint, 

s P — p 

— — = constant. 

p pT 

Since P = p, this tells us that the specific entropy of an isolated ideal gas is uniform. We 
have seen that this additional constraint is needed to identify the thermodynamic pressure 
—dF/dV with the gas pressure that appears in the Bernoulli equation. 

The function that we referred to as the hamiltonian actually deserves the name. In 
order to derive hamiltonian equations of motion we need to identify canonically conjugate 
momenta, in particular, a variable that is conjugate to the variable p. As shown by Fetter 
and Walecka (1980), a velocity potential fills the role admirably, with v = — grad$ by 
definition. The continuity equation takes the form p = the Bernoulli equation 

(in integrated form) is $ = 5H/5p. 

II. 3. Dynamical action principle 

The ideal gas is a system governed by an action 

A = J dt J d 3 x£[p,T, <$>,k ,P], 

where E C 1R 3 is connected, ko G H, p, T and $ are scalar fields on E and P is a scalar 
field on the boundary <9E of E. 

For the present, the domain E is fixed and the external pressure P is irrelevant; the 
action is varied with respect to the density p, the temperature T and the velocity potential 
<E>. The irrotational velocity field v = — grad $ is subject to the boundary condition 

v\di; is tangential. (2-17) 

The lagrangian is 

C = p$ - h = p(§ - v 2 /2 - (/) + X) - V(p, T, S), 
where v = —grad $ and A is a lagrange multiplier related to the fixing of the mass, 




The potential V is 

V(p,T,S) = f(p,T) + sT, S = s/p= -TZlnko = constant. 
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Variation of T gives Eq.(2.13), 



dV/3T = 0, ^ + s = 0. 



, P 

In — ; = n 



For the ideal gas it is 

In 

T n k 

Variation of $ gives the equation of continuity (local conservation law) 

p + div(pv) = 0. 

With the boundary condition (2.17) imposed on v, it ensures the global conservation of 
the mass. Variation of p leads to 

$-v 2 /2-cP + \= ^. 

dp 

Taking the gradient we get, since the entropy is uniform, the differential form of Bernoulli's 
equation (Bernoulli 1738), 

p— + pgrad(p = -grad p, 

with the gas pressure 

- V. (2.18) 



dV 

v ■= p^r- 

op 



T,S 



Besides the mass, the only globally conserved quantity is the hamiltonian H = J E d 3 x h, 
with the density 

h = p(v 2 /2 + (j)) + V[p]. 

The local conservation law is 

h + div(v(h+p)) = 0. 

The parameter ko is free, independent of volume and mass. That a free parameter must 
appear is evident since it must be possible to change the configuration of the system by 
heat transfer. 

The lagrangian, for each fixed choice of ko, desribes configurations of an ideal gas 
restricted to an adiabat. The equations of motion preserve mass and energy and the 
system is thus in a sense isolated. It is already possible to discuss the "changes" that form 
the central subject of thermodynamics. 



II. 4. Cooling 

By "cooling" is meant a slow loss of energy, as the system goes through the equilibria of 
a sequence of adiabatic lagrangians indexed by ko, this being the only parameter available. 
It is regarded as a reversible process, for it is supposed that heat can be supplied, slowly 
so as to make the departure from equilibrium negligible, to reverse the cooling process. 
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The domain E is supposed fixed, p and T uniform and $ = 0; hence p is fixed, unchanged 
during the process of cooling. 

During this process, the on shell internal energy density is u = nTZTp, the total energy 
is U = nlZTM, and changes in U, T and ko are related by 

dU = nTZMdT, dk = d-£- = -nk ^. 

As this represents a heat loss, In ko is increasing. Hence 

dU = -TZTMdlnko := TdS, 

with 

S := -TZMlnk . (2.19) 
This relates the parameter ko to entropy. We define the entropy density by 

s = —TZplnko- 

The specific entropy density s/p is uniform in this case. From the point of view of thermo- 
dynamics, with its primary emphasis on equilibria, the lagrangian density is a function of 
the 3 independent variables p,T and feo- The on shell condition (2.5) is the fundamental 
relation of the ideal gas: 

s = -Tip (ln^-n). 

A dynamical description of cooling as an ongoing process may rely on the heat equation. 

II. 5. Free expansion 

It is a process that begins (at t = 0) with a configuration in which the gas is uniformly 
distributed in a fraction a of the volume, with density p and temperature T, at rest. From 
these initial conditions the gas spreads, eventually occupying the total volume. Thermo- 
dynamics deals mostly with equilibria, but the adiabatic lagrangian allows for motion as 
well, and it describes the agitated configurations of the expanding gas for t > 0. Since the 
equations of motion preserve energy, this predicted, adiabatic motion of the gas cannot 
lead to equilibria. To bring the gas to rest another process is required, involving a change 
in the value of ko- 

It is observed that the gas comes to rest, and for this to happen a non adiabatic process 
must take place, with a change in entropy. In order that the time development described 
by the adiabatic lagrangian retain its relevance, one must postulate that the additional 
process take place on a much longer time scale. 

We thus envisage an initial, rapid, adiabatic expansion, leading to configurations of 
non zero kinetic energy, and a slow conversion of that kinetic energy into heat. Without 
attempting to describe the latter, it is nevertheless possible to predict the end result. If 
the energy of the asymptotic, equilibrium configuration is the same as that of the initial 
configuration, then since both are equilibria of the same ideal gas we have U = nlZTM = 
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U' = nlZT'M, the prime referring to the asymptotic state. Hence T = T' , and since 
p' = ap, k' = ako and the change in the function S is 5S = — TZMlna, which represents 
an increase since the fraction a is less than 1. In general, if some heat is lost or supplied, 
we shall have 

5S := -TIM 5 In k = TIM In k ° 



and on shell, 



o V 

5S = KM(\n^ +nln— ). (2.20) 



II.6. Work 

We consider changes of volume, and begin with the case of a cylinder with a piston. 
The cross section is A and the longitudinal variable is z, with < z < z\. Conjugate to 
z\ is the external pressure P excerted by the piston. 

A change dV = Adz\ of the volume requires that the gas supply energy in the form of 
work, dE = PdV = PAdz\. The total, conserved energy of the gas is thus modified by 
a term +PV (an additive constant has no significance) and the lagrangian is modified by 
the inclusion of a term — PV = —PAz\ , 

/ d 3 xC^ [ d 3 x(p(^> -v 2 /2 + X) -TZTp In— - PAz u 
Jt, Jt. v T n k J 

The action is to be extremized, in particular, by variations that vanish at <9£. The earlier 
equations of motion therefore continue to hold. In addition, we must consider the one 
parameter family of mass preserving variations of the form dp/p = —Szi/zi, 5$ = 5T = 0. 
This gives an additional equation of motion, 

/ d 2 xC -—I d s x (p($ -v 2 /2 + X)- TZTp{\n + 1)) - PA = 0. (2.21) 

The first term comes from changing the boundary, the second from variation of p and 
the third from the new term that was added to the lagrangian. With p = TZTp from the 
equations of motion; after multiplication by z\ it is 



z x [ d 2 xC- [ d 3 x(C-p) - PV = 0. 



By the equations of motion, C = p. The conclusion is that the external pressure P is equal 
to the average of the gas pressure over the face of the piston. 
The hamiltonian associated with the expanded lagrangian, 



H = J d 3 x (pv 2 /2 + UTp In JL-) + PV 
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is identified, in the case of an equilibrium configuration, with the Helmholz free energy 
U + PV. In general, the extra term that is needed in the case of a variable boundary can 
be expressed as 



where P is a field that coincides with the external pressure at the boundary <9E. 
In open atmospheres P = at the upper end; it implies that T = there. 

Remark. In Einstein's theory of gravitation the cosmological constant can be interpreted 
as an external pressure on space time at infinity. 
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III. Mixtures 

III.l. Basic postulates 

The adiabatic action for a system of two parts, occupying distinct portions of a total 
domain E, with no interaction or mutual constraint between them, is the sum of the 
adiabatic actions of each, 

/ C= I dd 3 z+ / £ 2 d 3 z, EiUE 2 = E. 
In the case of two ideal gases, 

j^ = Y,j^ ^(piC^i-^/a-^ + AO-v;), 

with 

V i = 1l i Tp i \n^, ki = ^, i = l,2.. 

The temperature field is defined over E; we may write Tj = T| , i = 1,2. 

We restrict ourselves, temporarily, to the approximation in which the gravitational 
field is ignored. In this case there is a unique, static solution of the equations of motion 
with uniform densities pi, P2 and temperatures Ti,T 2 given by pi = Mi/Vi and (eTj) ni = 
Pi/koi, i = l,2. 

We want to know what happens to the system if at t = the barriers that have confined 
the two gases to their respective domains are removed. Each gas will expand (diffuse) into 
the region originally occupied by the other. This implies motion and kinetic energy. Since 
energy is conserved, the adiabatic equations of motion do not lead to a state of equilibrium. 
But a real physical system can be depended upon to find its way to rest, it is therefore 
necessary to postulate some non adiabatic process, with a longer time scale. As in the case 
of simpler situations examined in Section II, and for the same reason, we may suppose 
that this slow process takes the system through a sequence of lagrangians indexed by k i 
and fco2, the only parameters that are available. The changes can be interpreted in terms 
of a changing (increasing) entropy. 

We expect, then, that the mixture will eventually be governed by the same lagrangian 
(but both densities will extend to the whole domain), with a new set k' Q1 , k' Q2 of parameters. 
Can we predict these new values? More precisely, the problem is to determine a relation 
between the two parameters. As heat is added to or withdrawn from the mixtures the 
entropy changes, tracing a path through the fcoi, fco2 plane. We need to fix this path. 
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Fig 1. Example of a path followed by a two component system in entropy space. The 
coordinates are the entropy parameters fcoi and &02 and the shown shape conforms to the 
case of disassociation of atomic hydrogen. 



The Gibbs-Dalton hypothesis (Gibbs 1876, Dalton 1808), justified within the kinetic 
interpretation of the ideal gas as an ensemble of non interacting particles, predicts that 
each gas undergoes the same change of specific entropy as would be the case if the other 
gas were absent. It is assumed that no chemical reaction is taking place. 

The mathematical statement of the Gibbs-Dalton hypothesis is, in view of Eq. (2.20) 

h' n' T- 

In 221 = l n ^ + ni l n il 7 

hi Pi T' 

where pi and are the densities and temperatures before mixing, p\ and T' are the 
densities and the common temperature of the final equilibrium configuration . 
The pre-mix parameters satisfy the relations 

Pi 



In — = m, i = 1, 2. 



It may happen that T\ = T2, and we may wish to equalize the temperatures before mixing. 
But in any case the Gibbs-Dalton hypothesis is the condition that the two equations 



ln-^— =n l? i = l,2, (3.1) 

with equilibrium values of the densities, are solved by the same value of T' . A possi- 
ble interpretation is that the difference between the two temperatures defined by (3.1), 
dissipates with adjustment of the entropies. 

The extension of this constraint to non uniform configurations is problematic, as will 
be seen. 
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This Gibbs-Dalton constraint imposes a relation between k' Q1 and k' 02 , 

(k' 01 /M 1 r = (k' 02 /M 2 r. (3.2) 

To the extent that the Gibbs-Dalton hypothesis is verified experimentally we can regard 
this constraint, and the lagrangian determined by it, as being dictated by experiments. 
The lagrangian of the mixture is now determined by a single, real parameter that serves 
as a record of heat that is lost or supplied. 
A natural choice for the lagrangian is 

C = J2 (Pi@i -v t 2 /2-<p + Xi) - TZiTpi In ^-j , (3.3) 

with pi, p 2 defined over E, the masses Mj = J E d 3 xpi fixed, and constants ha satisfying 

(W^i) n2 = (W^2) ni - (3.4) 

We have dropped the primes on these parameters. 

The equations of motion are the equations of continuity, the Bernoulli equation, namely 

$ l -v l 2 /2-<f> + \ l =n i T(ln^ + l), h:=^, i = 1, 2, 

and from variation of T the adiabatic condition 

Y / K i p i (\n^-n i ) = Q. (3.5) 

Consider the uniform distribution, for which pi = Mi/V. In that case, because of (3.4), 
both expressions In ki/koi — rii, i = 1,2, vanish for one and the same value T ' of T. We have 
Inki/koi — rii = rii ln(T ' /T) so that the third equation reduces to RiPi n i ln(T'/T) = 0, 
with the unique solution T = T' . Hence Inki/ka — rii = 0, % = 1,2, at equilibrium and 
both gases are polytropic. 

Before mixing, both terms in (3.5) vanish separately on shell. The constraint (3.4) 
ensures that there is a solution for which this is true for the mixture as well. 

The total pressure; that is, the pressure on the walls, is 

dC 

P = C-J2prj r =Tj2^^- (3-6) 

OPi 

Individual pressures are suggested by inspection of the "Bernoulli equations" of motion, 

grad (vi 2 /2 + (f)+ — ) =: grad (v t 2 /2 + </>) + — gradp, = 0. 
opi Pi 

This yields pi = TZiPiT for equilibrium configurations, but nothing that is useful in general. 
Additivity of pressures is valid in the sense of (3.6), where each term is the pressure of a 
single gas oppupying the total volume. It is also valid, but only at equilibrium, in terms 
of the pressures that are identified through their role in the Bernoulli equations. 
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III. 2. About the Gibbs-Dalton hypothesis 

The constraint (3.4) is an expression of the Gibbs-Dalton hypothesis. This hypothesis 
is widely used in the context of equilibrium thermodynamics. It is, however, nonlocal. 
Consider a long tube filled with a mixture of two ideal gases. One can wait until mixing is 
complete and the temperature has become uniform over long distances. But if something 
should perturb the system locally, even if the disturbance is uniform over a moderate 
interval, then one is at a loss to understand how the constraint, depending as it does on 
the total masses, is to be implemented locally. We feel that a localized version of the 
Gibbs-Dalton hypothesis may be appropriate in such cases. 

A local version of the constraint would replace the masses by the local densities; re- 
quiring that, at equilibrium, both terms in Eq.(3.5) vanish separately. But this would be 
inconsistent with the other equations of motion whenever external forces, such as gravity, 
are present. 

The strongest version of the Gibbs-Dalton hypothesis would relate the two densities to 
each other under all conditions, not just in equilibrium configurations, thereby reducing 
the number of degrees of freedom. It seems that this is done routinely in some applications, 
which explains the fact that dynamical processes have been described in terms of a single, 
uniformly mixed density distribution. (Exception: superfluid helium.) 

The truth may lie somewhere between these extreme interpretations of Gibbs-Dalton; 
between a local version that admits only one independent density degree of freedom, and 
the global version with two densities that are locally independent while only the masses 
are involved in the constraint. Implications for equilibrium configurations are identical in 
all versions, in the absence of external forces. 

One way to think about this problem may be to postulate a slow, non adiabatic process 
that is observed as a tendency of the mixture to anull the difference pi/fcoi — /°2Ao2- 
It is an irreversible process that calls for a role to be played by entropy. A Lagrange 
multiplier may be interpreted as the limiting case of a force that implements this tendency 
instantaneously. But it is characteristic of all irreversible processes that they are slow, and 
they can be taken into account in dynamical processes only under that condition. The 
local version, with its strong identification of densities, is difficult to justify. Only global 
versions will be studied in this section. 

III. 3. Dissociation 

We consider the simplest example of a chemical reaction involving ideal gases, the 
transformation between atomic and molecular hydrogen, 

H 2 <-> 2H. 

We fix, once and for all, the domain E with volume V and the total mass M. 

At sufficiently low temperatures the gas is almost purely molecular, with molecular 
weight and adiabatic index 

fii = 2, n± = 5/2. 
At equilibrium we shall have v = and 

Pl = p = M/V, In— — = m. 
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At sufficiently high temperatures the gas is almost purely atomic, with p 2 = 1, n 2 = 3/2. 
At equilibrium, 

p 2 = p = M/V, ln^- = n 2 . 

Consider a mixture of two gases, similar in all respects to H and H 2 , except that the 
reaction analogous to H 2 -B- 2H\ does not take place. If we assume the Gibbs-Dalton law 
for this mixture then we are led to the lagrangian density 

£ = J2 (f>i(*i - ^V2 - (f> + K) - KiTpi In 7f ^-) , 

with values of the parameters /coi and /C02 that ensure the compatibility of the on shell 
conditions \n(pi/T ni k i) = rii, i = 1,2 that must hold at equilibrium. A hypothesis that 
can be tried is that this Gibbs-Dalton lagrangian retain some validity for the real hydrogen 
problem. But it must be changed in essential ways. 

In the case of real hydrogen the individual masses are not conserved; there is only one 
equation of continuity and only one velocity field. Consider therefore 

J £ = J^d 3 x (p(* - v 2 /2 -<p + \) + e Pl -J2 KiT Pi In ^^-) , 

with p = pi + p 2 and e > constant. This expression, with constant parameters ka, is a 
reasonable candidate for hydrogen, but the Gibbs-Dalton constraint on the parameters is 
no longer applicable. 

The term ep± is the binding energy of molecular hydrogen. Variation with respect to T 
yields the relation 

E^( ln T4o~-^)=0 5 (3-8) 

a relation that we shall put aside for now. 

Variation of the densities with p = p\ + p 2 fixed, in the case of equilibrium, gives 

d{^ n ^ Tln T^k^~ €pl ^ = °' 

or explicitly 

ftiT(ln— ^— + 1) -TZ 2 T( In— ^— + l) = e. (3.9) 

This relation appears in textbooks on thermodynamics, except that, in some of the books 
consulted (Holman 1974, DeHoff 1993), the binding energy is not taken into account. And 
the unknown entropy of a reference configuration used in these books is here represented by 
the term IZiTpi In koi in the lagrangian, parameterized by koi, fco2- Thes parameters, and 
thus the entropy of the reference state not unknowable; but determined by the experiment. 

Relations (3.8) and (3.9) are the equations of motion that determine the configurations 
of adiabatic equilibrium, for each choice of the parameters koi and ko 2 , unique equilibrium 
values of T and for the ratio pi/p 2 . (The total density p = pi + p 2 and the volume are 
both fixed.) 
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Addition or withdrawal of heat, or a change in energy by other means, produces a non 
adiabatic path through the 2-dimensional space of lagrangians with coordinates &oi, ko 2 . 
To discover the relation between T and pi/p 2 we must know this path. See Fig. 1. 

Since 'R, 2 = 27^i, Eq.(3.9) reads 

ln(^) - (m - 2n 2 ) InT = e/K x T + 1, q := ^1 , 



P2 



hi 



or 



_^ T l/2 e - e /R!T. 



(3.10) 



ppi 1 — r ep 

If the parameter g is constant, then this is a special case of the famous Saha formula 
(Eggert 1919, Saha 1920). That q is constant is thus another property of ideal gases that 
is suggested by statistical mechanics. It is also supported by observation. The molecular 
binding energy is about 4.5 ev, or in relation to the rest energy, in energy units e = 
(4.5/1876 x 10 6 )c 2 « 2.15 x 10 12 . Thus 

e/Kx = 4.3 x 10 12 /. 83214 x 10 8 = 51674. 

Fig. 2 shows r versus T, for ni = 5/2, n 2 = 3/2 with q = 1. This result shows that a 
reasonable lagrangian for hydrogen across the dissociation region is obtained by assigning 
a constant fixed value to q. 
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Fig. 2. The ratio r against T, Eq.(3.10). 

The measurements are usually performed under conditions of constant pressure. Setting 
pi/2 + p2 = p/lZT we get, when n 2 = 5/2 and n 2 = 3/2, 



^-P = T^~ 2 P = ^ T 3 / 2 e-/^ T ; r = p 2 /p. 
ppi 1 — r z 2e 



(3.11) 
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This relationship is plotted in Fig. 3 for qlZ/p = 1. 
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Fig. 3. The ratio r against T, Eq. (3.11). 

Eq.(3.10) applies to laboratory conditions of constant volume. Eq.(3.11) applies under 
similar conditions when the gas is kept under constant atmospheric pressure. Both equa- 
tions can be used for atmospheres if the dependence of temperature and density (resp. 
temperature and pressure) on altitude are known. 

Eq. (3.8) gives us other kinds of information. This equation allows to calculate both 
koi and /C02 in terms of the equilibrium values of pi,/92 and T. With (3.10) it allows to 
calculate all the thermodynamic potentials as functions of the same variables. For example, 
the expression for the on shell internal energy is the same as for the case of a mixture, 
corrected only by the addition of — ep\. 



III. 4. Other reactions 

We consider, briefly, a chemical reaction involving 4 gases, by which a uniform mixture 
of gases of types 1 and 2 are transformed into gases of types 3 and 4, according to the 
chemical equation 

i=l,2 i=3,4 
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The coefficients are the ratios of reagents, in grams. Three masses are conserved, v 2 M\ — 
v\ M%, V3M2 + v 2 M 3 , V4M3 — Z/3M4, and the three corresponding densities satisfy conser- 
vation laws. The kinetic part of the lagrangian density therefore must take the form 

Adn = (>2Pl - flP2)($l - Vl 2 /2 + Ai) 

+ (>3P2 + ^2P3) (^2 - V 2 2 /2 + A 2 ) + (Z/4/O3 - ^3P4) ($3 - #3 2 /2 + A3) 

and 

C = Adn + P<f> + epi ~ Y] UiPiTln 1 , 

l<i<4 m 

where p = Y^ Pi an d the term ep\ represents the energy needed to make the reaction go. 

At equilibrium the densities are fixed by the volume and the masses, up to a variation 
of the form 

{dpi, dp±) oc (z/i, z/ 2 , -^3, -v±). 
This, and the variation of T, gives the two equations of motion 



and 

E rf /^( ln + !) = dpie/^iT. 

1<<4 1 °* 

The latter can be expressed as 

„^l/Ml ^2/M2 7 ^l/Ml I, ^/M2 

In Pl . p2 . = ln 9 + (m + na - n 3 - n 4 ) InT - e/T^T, q = fco1 , 02 . 

_^3/M3 ^^4/^4 7 i,^3/M3 7,^4/M4 

P3 P4 ^03 ^04 

one recognizes the law of mass action (Guldberg and Waage 1864, Gibbs 1875, Holman 
1974, page 489), except that here the right hand side is not an unknown function, but a 
function that is known explicitly in terms of the 4 entropy parameters fcoi, &04. This 
formula is formally the same as one that one derives from kinetic theory, it provides us 
with a hint that the parameter q may be an invariant. 

It is tempting to speculate that q is a constant, in that case there remain 3 independent 
entropy parameters. It seems likely that they are subject to constraints of the Gibbs-Dalton 
type. Thus 

(W^i)™ 2 = {k 2/M 2 ) n \ (k 03 /M 3 ) n * = (k 04 /M 4 ) n *. 

In that case the index of entropy is just one free parameter, as it should be for a complete 
determination of the entropy. 
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IV. The propagation of sound 



IV. 1. The ideal gas 

The dynamical equations that govern the adiabatic excitations of an ideal gas along a 
fixed direction are (Section II. 3 and Laplace 1816) 

p + ( pv y = o, ^-vy 2 + x = nT(\n T ^- + i) (4.i) 

and 

l n A =n fc :=_£_. (4.2) 

For perturbations around static equilibrium, to first order, we have 

k 

dp = -v'p, dk = 0, v = -n(ln — + l)dT'; . (4.3) 

ko 

Thus 

dp/p = -v' = K(n + l)dT", dT'/T = ndp'/p, 

and finally 

dp/p = HT(1 + -)dp"/p. 

n 

The speed of sound is thus y/IZT'j, 7 := 1 + 1/n. Only one degree of freedom is excited. 

Consider next the implications of the global version of the Gibbs-Dalton constraint. In 
a mixture of two ideal gases two degrees of freedom are excited and two normal modes of 
propagation will appear. 



IV. 2. Propagation of sound in Gibbs-Dalton mixtures 

The linearized equations for a perturbation of a Gibbs-Dalton equilibrium configuration 
are 

r)V 

*'=-" :=i %U (4 - 4) 

dT" __ R^p'j + R 2 dp% 

t ~ ni n lPl + n 2 n 2 p 2 ' ^ ' ' 

The second equation, - from Eq.(3.5) - defines an adiabatic change. It is assumed that 
the unperturbed densities and the unperturbed temperature are uniform and that the 
perturbations depend on only one of the spatial coordinates; the prime denotes the spatial 
derivative. 
Let 

k = Tl 2 /TZ 1 = px/ p 2 ; r = pi/p 2 . 
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then 



(4.6) 



Next 



dm TIT t dp" TIT k dp'{ TIT dp'l 

— = J— + J— + — • 4 - 7 ) 

pi pi A pi pi A p 2 pi pi 

dp 2 _ 1ZT k dp'{ + KT L dpl + TZTdpl (4 g) 



p 2 P2 A pi p 2 A p x p 2 p 2 

The matrix of squared velocity is 

TIT ( 1 + 5 f 

A 2 , A = mT + n 2 K. 



For mixtures we define the reduced speed c in terms of the cgs speed v = cy/lZiT. This 
reduced speed is a zero of the determinant of the matrix 

1 + - - c 2 - \ 

I $ 2 J ( 4 - 9 ) 

The theory predicts two modes, with speeds that vary with the concentrations. Experi- 
menters report a single mode. We shall see that this discreapancy can be overcome, but it 
is of some interest to begin by describing the predictions of the naive, Gibbs-Dalton model 
in some detail. 

Two different types of binary mixtures must be described separately. 

Type 1 mixture, "similar" gases 

This is the case when K72 > 1, for example: 

Nitrogen/ Argon: k = 28/40, m = 5/2, n 2 = 3/2, 

with K72 = 35/24. There is a Laplace mode, with speed intermediate between the adiabatic 
speeds of the two pure gases. The two amplitudes are in phase. The dominant gas carries 
most of the energy but the amplitudes are comparable. A second mode has a speed that 
interpolates between the Newton (isothermal) speeds of the two pure gases. The two 
amplitudes are in opposite phase. The amplitude of the dominant gas tends to xero in the 
limit when this gas is alone. 

All this tends to be confirmed by experiment. The second mode seems to be anti 
intuitive and there are several reasons why it may be expected to be unobservable, or at 
least to justify the fact that it has escaped detection. The prediction of the Gibbs-Dalton 
model needs only a minor correction to agree with experiment. See Fig. 4. 
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Fig 4. Speed of the two modes of sound in Ni/Ar, plotted against hi^-. The curve just 
below the data points is the speed given by Eq.(4.10). 

Type 2 mixture, "disparate" gases 

It is the case when K72 < 1, for example, 

Helium/ Argone, k = 4/40, ni = = 3/2. 
There is a mode, with speed that approaches the adiabatic speed of sound in the lighter 
gas in the limit when this gas is alone (NE corner), but in the opposite limit it approaches 
the Newton value of the heavy gas (NW). The amplitudes are in phase. 

A second mode has a speed that approaches the adiabatic speed of sound in the heavier 
gas in the limit when this gas is alone (SW). Surprisingly, the amplitudes are in opposite 
phase; the two gases are moving in opposite directions. See Fig. 5a. 
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Fig. 5a. The speed of sound in He/ Xe, a type 2 (disparate) mixture. Coordinates as in 
Fig. 4. The interpolating curve is from Eq.(4.10). Data points from dela Mora (1986) and 
Bowler (. 
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Fig. 5b. The speed of sound in He/Ar, a type 2 (disparate) mixture. Coordinates 
as in Fig. 4. The interpolating curve is from Eq.(4.10). The high datum point is from 
Smorenburg 1996. 

The physical interpretation is more difficult for this type of mixture. Experiment yields 
a single mode with a speed that varies with concentration and that approaches the expected 
values in the limits of either pure gas. It appears that a cross over takes place; in fact, 
there is evidence of a "critical" concentration. 

Experimenters report a single mode in all cases, with a speed that varies smoothly with 
the concentration. There are frequent hints, however, that there are other, less prominent 
modes, about which no information is given. 

The reason for this disagreement was at first attributed to the fact that no account had 
been taken of any kind of damping, and for high frequencies this has been confirmed. But 
at low frequencies the loss of energy to absorption is not thought to be important. Indeed, 
the measured rates of absorption are extremely low. (Holmes 1960) 

The formula that fits the observations, and much better than could be expected, is 



The origin of this prediction is not thermodynamics but calculations of particle dynamics 
and the Boltzmann scattering equations. For a brief account of the calculation see de la 
Mora and Puri (1986). It is predicted to hold in the limit when the effect of diffusion is so 
strong as to force the two velocities to be nearly equal. That is the key to the success of 
the formula, and that must be the basic fact that is responsible for its success: apparently 
something forces the two velocities to be nearly equal. But the conclusion that it is due 
to strong diffusive damping is not compelling. 

The equations that have been proposed to explain the success of (4.10) are equations 
for the two velocities and, some times, two temperatures. We have noted that the intro- 
duction of two temperatures may be interpreted in terms of a deviation from Gibbs-Dalton 
equilibrium conditions. But we are strongly constrained by the need to retain the standard 




(4.10) 
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equations, continuity and Bernoulli. Indeed, any modification of the kinetic part of the 
lagrangian will cause a lack of conservation of masses: that is, of particle numbers. It is 
therefore of interest to explore alternatives. 

The data quoted are from experiments at frequencies around 1 MHz, except for that of 
Smorenburg (...), the single, isolated point at a high speed in Fig. 5b, where the frequency 
was of the order of 10 9 Hz. More data is needed, over a wide range of frequency and 
temperature. 



IV. 3. A simple model of interactions 

Let us add the following "interaction term" to the potential; that is, to the hamiltonian, 

(4-H) 

with a constant. The reason for this choice will become clear almost immediately. Such a 
term will affect the formula for the internal energy, but it does not change the formula for 
the entropy. The hypothesis of Gibbs and Dalton can be maintained, and we shall assume 
values of the constants accordingly. The formula (3.6) remains valid; it is the Dalton 
rule of additivity of pressures. 

The equations of motion that would be recognized as Bernoulli equations, 

grad (*< 2 /2 + 4> + KiT{ In + l) + | 

no longer suggest a useful definition of individual pressures. The acceleration of the 
molecules of each gas is strongly dependent on the configuration of the other. 

Again we consider first order perturbations of the equilibrium configuration. The lin- 
earized equations (4.4) are modified: 

d Pt = -v t = d x [ — 

Eq.s (4.7) and (4.8) gain additional terms, 

Fp2(dp2 dp" 







What makes the choice (4.11) special is the fact that, for a monocromatic wave, these 
additions are proportional to the difference v\ — v-i of the velocities. Comparison with the 
standard theory shows that this is what is needed. 

The reduced speed of propagation is now obtained by setting to zero the determinant 

l + ^-/3r-^-c 2 4+/3t"5 \ a 



^+/3r5 K + ^-p T h -c 2 J ' 4ftiT 



P ■■= TFT^F- ( 4 - 12 ) 
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In the limit of large (3 the only eigenvector has dp\j p\ = dp2/ P2 and the eigenvalue c 2 is 
precisely as in (4.10). 

Numerical results for a type 1 mixture are in Fig. 6. It may be seen that a value (3 = .5 is 
sufficient to bring the theory into substantial agreement with the experient. Higher values 
of (3 improves the agreement for low frequencies but tends to spoil it when compared with 
the few data that are available for the highest frequencies. 



1.05- 




Fig.6. The effect of the interaction (4.11) on sound speed in N%/Ar. The middle curve 
is from (4.12) with for (3 = .5. 



This modification of the lagrangian will be successful if the observation of a single, 
dominant mode is limited to low and moderate temperatures. If instead it extends to high 
temperatures an alternative may become more attractive (see below) . The theory predicts 
no dispersion. 

The case of type 2 mixtures (Fig. 7.) is more dramatic, very small values of (3 are 
enough to eliminate the isothermal modes in the SE and NW corners of the figure. The 
former disappears and the latter drops down to join the Laplace branch at the lower left. 
Agreement with experiment requires a value of (3 at least equal to 1, near perfect agreement 
is gotten with f3 = 5, and larger values of (3 only improve the fit, except for the data at 
the highest frequencies. 
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Fig. 7. Effects of the simple interaction term on the speed of sound in He/Xe, Eq. (4.12). 
From top to bottom: f} = 0, (3 = .01, = .1,0 = 1; 



IV. 4. A variant 

One may try to try to replace a by aT. The potential is then 

V = KiPiT In + aT^pTp- 2 . 

The new term may be considered as an addition to the free energy, and the expression for 
the entropy density is affected; it is possible to retain the homogenous entropy predicted 
by Gibbs-Dalton theory. Thus we retain the former values of the parameters hoi. 
Variation of the temperature gives 

The individual gases are no longer polytropic, even at equilibrium. The expression for the 
internal energy is independent of a and the sum formula (2.6) for the total pressure also 
remains valid. 

Adiabatic changes obey 

— ^ Ripm = ^ Hid* ((In + 1 j + - rfw ^ — • 

Define 

AT <V 1 Pi AT Nl ~ N * 

{eT) n *k 0i 2 
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then 

The linearized equation of motion is 

pi ~ Vl ~ dx d Pl > 

with 

^ = ™ l(ln r^7 + 1)+ 2\/^ 

The adiabatic variation is 

d^- = dT( Rl + *)+ k iT ^ + ^rM( d -^ - 

opi Pi Pi 4 y pi p 2 pi 

= (fi 1 + £)(£ft^ + Ar(!^-*l)) 

l^KiPiUi^ pi pi p2 

+ n t— + —T fp^( dp2 - dpi ) 

1 pi 4 V Pi P2 Pi 
This time, c 2 is an eigenvalue of the matrix 

1 + Ai(pi + N/Ki) - fir' 1 / 2 A lK (p 2 - N/n 2 ) + pr' 1 ' 2 
A 2 (pi + N/Hi) + fir 1 ' 2 k + A 2K (p 2 - N/Ki) - fir 1 / 2 

where 

1 /„ N \ 

^UiPini^ pi' 

and (3 = (a/4TZ\). For large a the eigen vector has dp\j p\ = dp 2 j p 2 and a is eliminated 
by adding pi times the sum of the elements in the first row and p 2 times the sum of the 
elements in the second row, with the result 

A{r iP 2 + n 1 T Pl + A 2 n 2P 2 + tz 2 t P2 - v 2 ( Pl + P2 ) = o, 

which again agrees with (3.10). 

Both models can account for the experiments within a range of temperatures and fre- 
quencies. The only observation at a temperature other that room temperature that is 
known to us is that of Smorenburg et al (Smorenburg 1995). At a molar concentration of 
3:1 in He/Ar (Type 2), (pi/p 2 = .3) , at 160K and 370 bar, they observe the same sound 
speed as in pure helium. This may indicate that, of our two models, our second one is 
closer to observation. 
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V. Application to atmospheres 



V.l. Introducing gravitation 



Gravitation enters as a component of most, if not all applications of physics. There is a 
set and tested procedure for introducing the gravitational field into any dynamical context. 
In classsical, non relativistic physics it consists of adding the gravitational potential energy 
to the hamiltonian. For an ideal gas one adds 



where p is the mass density and <fi is the Newtonian potential. The adiabatic lagrangian 
becomes 



When the effect of gravity is taken into account in the dynamics we get a theory of 
atmospheres. In the case of terrestrial gravity <p = gz, where g is a constant and z is the 
elevation. At equilibrium k is constant and the equation 



predicts a constant temperature gradient (lapse rate). The effect has never, to our knowl- 
edge, been tested in the laboratory (but see Graeff 2009), but the constant lapse rate 
is a feature that is observed in the earth's atmosphere and in the internal structure of 
stars (Lane 1870, Ritter 1880, Eddington 1926). For the earthly troposphere this formula 
actually gives a value of the temperature gradient that is closed to observation. 

We have noted that this successful application of a theory originally constructed to 
account for laboratory experiments where gravitational effects are insignificant, modified 
in standard fashion to include gravity, does not attribute the observed temperature gradient 
to the radiation from the sun or to any external source other than gravity (Fronsdal 2010). 

And yet it is evident that, in the absence of the sun, the present state of the earthly 
atmosphere could not endure. The extinction of the sun would lead to a general cooling of 
the atmosphere. This cooling is not described by any process encompassed by the adiabatic 
lagrangian dynamics. Instead, an external agent enters the picture, the spontaneous emis- 
sion of infrared radiation. As this is a slow effect we can regard the cooling as a sequence 
of equilibria of adiabatic dynamics, a slow loss of entropy and an increase in the value 
of ko. (Section II. 4.) The reverse effect is also possible and if the terrestrial atmosphere 
is stable over long times then we must conclude that spontaneous cooling by emission is 
balanced by heating provided by the sun. 





The gravitational field appears in the equation of motion, 



§-v 2 /2-4> + \ = 1lT(\n— + 1), k = p/T n . 



X-gz = KT(n + 1) 
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V.2. Mixed atmospheres 

At one time it was believed that, in an atmosphere consisting of several components 
with different molecular weights, the lighter gas would float on top. This was surely based 
on observation, since cooking gas, entered at ground level, tends to remain there; however, 
it does so only for a short time. Dalton made the radical proposal that each gas behaves as 
if the other were absent. This was an overstatement and led to much misunderstanding and 
debate, some of it recorded in Dalton's book (Dalton 1806). In the case of non evanescent 
gases that have had time to settle, Dalton's prediction is fairly close to the truth, but 
he did not go far enough, for actually it is observed that the concentrations tend to be 
independent of elevation. In order to account for this, one needs to invoke the intervention 
of a mechanism, called convection for short, or "mixing due to fluid motions" , or diffusion. 

This takes us back, once more, to consideration of a dissipative process, which can be 
interesting, but if the problem that concerns us is the density and temperature profiles of 
the mixed atmosphere then it begs the question. If, as is often observed, an equilibrium is 
ultimately reached, then we are mainly interested in the end result, and less in describing 
the process that leads to it. The process is dissipative and entropy producing; its prin- 
cipal feature is that it is slow. It can be ignored in the study of the final equilibrium 
configuration, and the adiabatic perturbations of it. 

We have embraced the Gibbs-Dalton hypothesis and we have incorporated it into our 
lagrangian for the mixtures that are not under the influence of gravity. But in the presence 
of gravity there is an ambiguity. We shall describe an attempt to adopt the global version of 
the hypothesis to atmospheres. It is only moderately successful. Then we shall determine 
the lagrangian by direct appeal to experience, to get a perfect fit. 



V.3. The Gibbs-Dalton model 

Consider at first a vertical column of air consisting of nitrogen and oxygen and confined 
to the range < z < z\ of elevation. Because the Lagrange multipliers are freely adjustable 
we can fix z = at ground level. If gravity were not present we should have an equilibrium 
configuration with uniform densities and temperature and the equations of motion would 
give us 

p, = /co*(eT) ni , Ai =7£iT(ni + l), i = 1, 2, 

implying the Gibbs-Dalton constraint, at equilibrium (pi/fcoi) = (P2/&02) and 

thus (M 1 A 01 ) 1/ni = (M 2 A 02 ) 1/n2 

Having chosen the domain and the parameters (except the Lagrange multipliers), we 
add the gravitational field. The masses, presumed unchanged, are now related to the 
partial pressures at ground level, * 

gMi = HiPi(0)T(0). i = l,2. 

The relative concentration of the heavy gas, pil p\ = (^2/ fJ>i)(M2/Mi) is increased at the 

lower level. We take the total density at ground level to be .0012, thetemperature300K ', themassratio3.26 

* It is assumed that the column extends upwards until the pressure (and the temperature) 
becomes zero. 
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( Pl (0) = 3.26p 2 (0)) and p x = 28, p 2 = 32, Pl (0) = .9184 x 10" 3 , p 2 (0) = .2816 x 
10" 3 , m=n 2 = 2.5. 

The constraint (3.4) gives, for the masses implied by these parameter values, 

k 2/koi = (/xi//x 2 )(l/4) = .26825. (5.1) 

The equations that determine the equilibrium configurations are 

\ 1 -<f> = K 1 T(hix+l + n 1 ), x := ^ - , (5.2) 

(eT) ni fe i 

A 2 -0 = ft 2 T(lny+l + n 2 ), y := (5-3) 

Tli P i\nx + TZ 2 p2^y = 0. (5.4) 

It is useful to note that y(0)/x(0) = p, 2 / P-i- The third equation can be rearranged to give, 
since n\ = n 2 , 

lns(0)= , —. ln^, lny(O) = — 1 ln^. 

Thus In x(0) = -.02824, x = .9722, In y = .1053, y = 1.111 and 

fc 01 = ^(0)/(eT) n x = 4.975 x 10" 11 , k 02 = p 2 (0)/(eT) n y = 1.339 x 10" 11 (5.5) 

and finally Ai = 3.09645 x 10 9 , A 2 = 2.81343 x 10 9 . 

The atmosphere ends where the pressure vanishes; that is, where T = 0. Since Ai > A 2 , 
this happens at the point where 

= Ai, x = l, y = 0, z = Xx/g w 31.6 x 10 5 (32 km). 

Just below this point we find that dT oc dz, p\ oc T ni , p 2 oc T n2 e _Q / T , with dz = \\/g — z 
and a = (Ai — A 2 )/7?. 2 = 121.885. The first gas is dominant and adiabatic as if it alone 
constituted the atmosphere; the second gas has a different profile in this region. 

More interesting are the temperature and density profiles near the ground; we have 

T' 

- — (0) = 3.2303 x 10" 7 = 3.23 %/km. (5.6) 

In a nitrogen (oxygen) atmosphere the numbers are 3.14% (resp 3.59 %). The density 
profiles are 

Bl(o) = -7.852 x 10" 6 , ^(0) = -9.038 x 10" 6 . (5.7) 

Pi P2 

The last two ratios are found by observation to be very much closer. The concentrations 
are reported with four significant digits (no error bars) at all levels where the theory can 
be applied with any confidence (up to 10 km). 
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The slopes were calculated as follows. From (5.2-3), 

-gm = TZT^- + TZT'(\nx + 1), -gfi 2 = TIT— + TZT'(\ny + 1) 

Pi P2 

From (5.4) and some reduction, 

T_ __ -g\i\H2 piQnx + 1) + p 2 Qny + 1) 

T ~~ TIT ^ip 2 [(lny+l) 2 + n 2 ]+At2Pi[(lnx+l) 2 + ni]' 

The prediction (5.7), of two different rates of decrease of density, is not very far off. It 
may be correct when the rate of diffusion, that tends to reduce the difference, is small, and 
for young mixtures, when diffusion did not yet have enough time to act. The difference is 
only 1 percent per km, but it is much too large to be reconciled with measurements in our 
atmosphere. 



V.4. Entropy from experiment 

Let us take it as a result of observation that the concentration of the principal gases in 
our atmosphere are very nearly constant at low altitudes, and use this datum to find the 
values of the entropy parameters fcoi, 

The 3 equations of motion are 

y^RiPilnxi = 0, Xi := f P * , (5.8) 

\i- gz = Tl i T(lnx i + l + n i ), i = l,2. (5.9) 
We eliminate the Lagrange multipliers by differentiation, to get 

-g = H i T'(hix i + l) + H i Tp' i /p i , i = l,2, (5.10) 

and introduce the experimental evidence in the form p[/pi = p' 2 / 1 Pi- 

We begin the reduction of these equations by deriving 2 relations without logarithms. 
The first is obtained by differentiation of Eq. (5.8), 

p'i p' 2 T' Y^Kipirii 
- " — n— , n — 



Pi P2 T Yl Kipi 

The second relation is obtained by multiplying (5.10) by pi, summing, and combining with 
(5.8), and using the last result, 

-9(pi+P2) =T / ^^p l (n + l), 

or 

T' = -^-, -:= = 1/28.7179. 

K(n + 1) fj, Y,Pi 
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Since the two densities are proportional, the quantities n and /i are uniform, and so is the 
temperature lapse rate. 

We are now in position to evaluate the entropy. Eq.s (5.10) give us 

lnx t = (n+l)(^-l). i = l,2; 
p 

also uniform, (values -.1026 and .3827, x = .9024, y = 1.466.) Finally, 

Xi 

with values 5.359 x 10" 11 and 1.011 x 10 -11 , compare (5.5) and /C02/&01 = -1887, see (5.1). 

The entropies that are determined this way, using the observed constancy of concen- 
trations, are not very different from those calculated on the basis of the Gibbs-Dalton 
hypothesis. 

When heat is supplied or withdrawn, without loss or gain in total mass, we note that 

koi = ^r(eT)- n '-\ 

Xi l\,i 

where pi is the partial pressure. Evaluated at ground level, the first factor remains constant, 
so that the system moves on a path 

q ■= kol +1 /k^ +1 = constant. 

The entropy density, given by the on shell condition (4.8) as 

s = ^TZiPi ln/coi = ^TliPiilnpi - m IneT), 

is a function of the densities, the parameters and the temperature. The temperature lapse 
rate (uniform) is 9.67 K per km and the proportional rate at ground level — T'/T(0) is 
3.23 x 10~ 7 or 3.2 % per km. The fractional density gradient is p'(0)/p(0)) = 8.075% per 
km. 



V.5. Interaction model 

It is common to all models that failure of the Gibbs-Dalton hypothesis (about equi- 
librium entropies) to account for the observation of (a single) sound speed in mixtures 
is attributed to the effect of interactions between the atoms of the two species. But the 
properties of mixtures at equilibrium do seem to be accounted for by this hypothesis. It is 
therefore of interest to discover whether the observed atmospheric concentrations can be 
accounted for without affecting the entropies. A model of sound propagation examined in 
Section IV. 3 does have this property. 

This first model has the appealing property of preserving the simplicity of the Gibbs- 
Dalton equilibrium; the equation of motion that comes from variation of T has a unique 
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solution such that the two terms in Eq.(4.4) vanish individually; thus x = y = 1. The 
interaction term (4.11) affects the density- variation equations. Instead of (5.2-3), 

M Pi 



A 2 - 



K 2 T{l + n 2 ) + 



a pi 
2 V 92 



This is incompatible with constant concentration if \i\ ^ p 2 . This model cannot, by itself, 
do the job, but of course it may work in conjunction with an ajustment of the entropy. 
The second model leads to 

X 1 -<t> = K 1 T(l + n 1 ) + ^M, 

2 V Pi 



A 2 - = U 2 T(1 + n 2 ) + 



aT fp~i 



and this is compatible with constant concentration for one value of a, namely j3 = a/472. i 
.004. 



5.6. The hydrogen atmosphere 

Saha (1821) used his formula in studies of dissociation in stellar atmosphere, with pres- 
sure and temperature data obtained from other sources, with results that agreed well 
enough with observation. But there was no study of these profiles within the context of 
the theory. Here, with the advantage of a complete formulation of the dynamics, we can 
calculate the profiles directly. 

To describe a hydrogen atmosphere we include the gravitational potential in the la- 
grangian. An isolated atmosphere is characterized by fixed entropy. The equilibrium 
relations (3.10-11) remain valid; they are the two equations that remain when the gravi- 
tational potential is eliminated from the equations of motion. If instead we eliminate the 
logarithms we get the simple result given on Eq.(5.11) below. 

What is needed in addition is a relationship between r and T for constant entropy; that 
is, for constant values of the parameters koi and ko 2 . The other equation of motion, Eq. 
(3.8), is the fundamental relation of the gas, for which there is no counterpart in Saha's 
approach. From it we can extract a formula for the density. But the most revealing result 
is obtained by eliminating the density p = pi + p 2 from the three equations of motion. The 
result is the following two simple equations of motion, 

X- gz + er = T < R(n + 1) >, (5.11) 

and 

(m - n 2 ) InT = - _ + ^ > + (1 + r)-^= + constant. (5.12) 



33 



where 



< R(n+1) >:= ^Riirii + l) — . 



Eq.(3.10) - resp. (3.11) - is a relation betweeen r and T that holds through a sequence of 
equilibrium configurations with different energies (assuming that q is constant), at constant 
volume - resp. at constant pressure. Eq. (5. 11-12), on the other hand, hold throughout the 
atmosphere, with fixed entropy. 

Eq. (5.11) is a smooth interpolation between the atmosphere of molecular hydrogen at 
low temperatures and the atomic atmosphere of high temperatures. There is an important 
change in the lapse rate across the transition region. 

The Lagrange multiplier A controls the total mass and the constant in the second equa- 
tion is the entropy; both are free parameters. 

From (5.11) and (5.12) on easily obtains profiles of temperature and densities. Further 
work in this direction is deferred. 



34 



VI. Conclusions 

Let us assess what has been done and, especially, what has not been done here. 

(a) It was, of course, well known that irrotational hydrodynamics can be formulated as 
an eulerian field theory with an action principl. Conservation of energy is a cornerstone 
of thermodynamics and all adjacent sciences. What appears to be new is including the 
temperature in the set of independent dynamical variables, with respect to which the action 
is an extremum. By this means all required information about a system is stored, once 
and for all, in the lagrangian. Traditional thermodynamics is an on shell projection; more 
precisely, a partial projection on the solutions of one or more of the variational equations. 

(b) Given the gas law and the expression for the internal energy of an ideal gas it is 
commonplace to calculate the entropy as a function of, say, density and temperature. This 
is a fundamental relation of the ideal gas, from which all its properties can be deduced 
(Callen 1960). What may be new is that there is an off shell level of theory at which 
S, T and p are independent variables, and that the fundamental relation is the on shell 
condition obtained by variation of the temperature; it is one of the equations of motion. 

(c) The equations that constitute the Euler Lagrange equations do not contain any- 
thing new, as far as the one component ideal gas is concerned. What may be new is the 
total reliance on the chosen lagrangian for all subsequent applications. Any expansion or 
generalization of the system or of the context must be done (so we say) by making what- 
ever changes that are necessary in the lagrangian, in a manner so as not to upset any of 
the applications that have already been successful. This last is a very strong restraint on 
invention. It is the main reason for pursuing this line of inquiry. . 

(d) Mixtures are usually studied by means of the Gibbs function, but such studies are 
always, as far as we know, hampered by the fact that this function is known only relatively 
to that of a reference configuration. It seems that the lagrangian gives more information 
than what has been extracted from the Gibbs function. We have taken literally, perhaps 
more so than is traditional, the additivity of lagrangians for composed systems. We have 
been led to a suggestion for the lagrangians for mixtures that are precise except for the 
values of a small set of parameters, in number equal to the number of components. This 
has greatly facilitated the task of using experimental information to pin down the entropy 
of mixtures, something that we have found to be unexpectedly easy. This is our best result. 

(e) The suggestion of "second sound" in a gaseous mixture is not confirmed by observa- 
tion. We have suggested that this signals a limitation of the Gibbs-Dalton hypothesis to 
serve as a guide for dynamics. It was shown that some of the observations can be explained 
in therms of an interaction of densities, without the assistance of kinetic theory. Further 
exploration of this model must wait for more data, over a wider range of frequencies and 
temperatures. 

(f) There are seeds of a controversy in all this. There is no doubt that all the results 
can be obtained without using an action principle. But, as was stressed already, the 
action principle is a powerful guide to further applications. Indeed, the action principle, 
combined with what is known about gravitation, may be at variance with the common 
conception of an isothermal equilibrium of an isolated gas in the presence of a gravitational 
field. For this reason we have continued the investigation, begun in the first paper of this 
series, of the classical arguments that seem to have convinced all our thermodynamicists. 
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Results have been relegated to the Appendix, to avoid giving the impression that this is 
all we are trying to say in this paper. We suggest that measurements of the temperature 
gradient in a supercentrifuge should be carried out. Such measurements should relatively 
simple; it would settle once and for all (via the principle of equivalence) the question of 
the isolated atmosphere. A profound reassessment would be required, whatever the result, 
for an isolated, isothermal atmosphere does not have a natural place in thermodynamics 
(Fronsdal 2010), while the existence of a temperature gradient is widely believed to be a 
violation of the second law of thermodynamics. 
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Appendix. The controversial "isolated" atmosphere 



If a vertical column of an ideal gas is isolated from radiation, both emitted and absorbed, 
then it may be reasonable to believe that both may be left out of consideration. That is, 
an isolated column of gas in the gravitational field would obey the equations of motion 
and exhibit a temperature gradient. This conclusion was reached long ago by Loschmidt 
(1876) and others. It was refuted by Maxwell(1868) and Boltzmann(1909), as follows. 

If we would grant that an isolated, vertical column of air in equilibrium were colder at 
the top and warmer at the bottom, then a paradox would arise. Consider this arrangement: 
A heat bath with the lowest point at z = is maintained with temperature T = T (at 
z = if not everywhere). Below it is suspended a vertical cylinder filled with an ideal 
gas, in thermal contact at z = with the heat bath, but otherwise isolated. The theory 
predicts that, while the gas at the top of the cylinder has temperature T , the bottom will 
be warmer, with temperature T\ > T . It is proposed to utilize this temperature difference 
to run an engine, taking heat from the bottom of the cylinder and returning it to the bath. 
The argument is circular, for a complete analysis would require a theory of heat engines 
operating between different levels of the gravitational field. Maxwell circumvented this 
objection by replacing the heat engine by a second cylinder, filled by another ideal gas, 
with different characteristics (molecular weight, density, adiabatic index). This second 
cylinder, at equilibrium, would have the same temperature T at the top, but a different 
temperature Ti ^ T\ at the bottom. Now provide thermal contact between the bottoms 
of the two cylinders and suppose that To < T\ < at t = 0. Then both T\ and T2 will 
begin to change, in the direction of equalization. If "heat flow" is defined in terms of the 
gradient of the temperature then, according to Maxwell, a permanent, closed heat flow 
will be established; apparently, a perpetuum mobile of the second kind. 

We believe that the force of the conclusion is strongly affected if it is shown that no 
energy flows. 

The perpetuum mobile of Maxwell is one in which no work is being accomplished, "noth- 
ing really happens", since the energy flow of the final, asymptotic configuration can be 
assumed to vanish without contradiction. But what is this final state? Does it in fact 
exist? We believe that, at late times, a stationary state will be reached, and we try to 
calculate it. 

To remove a slight complication, let us take away the heat bath but retain the ther- 
mal contact between the two cylinders at both ends, isolating the system. The original, 
adiabatic equations of motion cannot apply as they stand, since they require two different 
lapse rates. Another process is needed, similar to "cooling" , in which there is a change in 
entropy. To account for this additional process we use the model lagrangian 

C tot = C[^,p,T,...} + jC[k,a], 

with the gravitational potential included as before in the adiabatic part. The field a is 
canonically conjugate to T. The domain is the union of the two cylinders. The equations 
of motion include the equation of continuity, unchanged, and, at equilibrium, 

-0 + A i = 7£ i T(ln^ + l), h:= Pl /T n % i = l,2, 
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1Zpi(\n-^- - rii) = cAai, 
cAT = 0. 

It is important to emphasize that the entropy is not driven directly by gravity - that would 
imply an unorthodox gravitational interaction and probably a violation of the equivalence 
principle. When the thermal contact is made, between the lower ends of the two tubes, the 
temperature at that point will eventually settle at a value T intermediary between Ti and 
T 2 , resulting in a temperature gradient different from that induced by gravity in either of 
the two separated tubes. Let 

5T[ = {f-T t )/t, i = l,2, 

where £ is the length of the tubes, or more precisely the difference in elevation between 
the upper and lower ends. This increment drives the entropy, thus 

Kiim + l)pi xrr , 
c°i = y oT t . 

Because the process is reversible, a'" + a' 2 n must be zero, which leads to 

T' ~ —g- Pl + P2 



n lPl (n 1 + i) + n 2 p2(n 2 + iy 

where, to this order, one should interpret the densities as averages. This represents a 
reasonable interpolation between the two extreme cases in which one or the other cylinder 
is absent. 

There is no energy flow within the tubes and, since the temperature is continuous, no 
energy passes the boundary. 

If instead one would justify the belief that a heat engine can be made to run on the 
temperature difference between the lower ends of the two cylinders (Graeff 2007), then a 
similar analysis involving entropy needs to be made. 
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